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We discuss recent results on the evolution of unpolarized parton densities and structure functions in massless 
perturbative QCD. Present partial results on the next-to-next-to-leading order (NNLO) evolution kernels prove 
sufficient for reliable calculations at not too small values of the Bjorken variable, x > 10 -3 . One order more can 
be taken into account at x > 0.2. Inclusion of these terms considerably reduces the main theoretical uncertainties 
of determinations of a s (to about 1% at the Z-mass) and the parton densities from structure functions. 



1. Introduction 

Structure functions in deep-inelastic scatter- 
ing (DIS) are among the quantities best suited 
for measuring the strong coupling constant a s . 
They also form the backbone of our knowledge of 
the proton's parton densities, indispensable for 
analyses of hard scattering processes at proton- 
(anti-)proton colliders like Tevatron and the 
LHC. During the past two decades DIS experi- 
ments have proceeded towards high accuracy and 
a greatly extended kinematic coverage [Q . 

To make full use of these results requires 
transcending the standard next-to-leading order 
(NLO) formalism Indeed besides the QCD 
/3-functions to even NNNLO |,[§, the NNLO 
(2-loop) coefficient functions for DIS have been 
calculated some time ago [||6|. However, only 
partial results have been obtained so far for the 
corresponding 3-loop splitting functions [Q— 13 ] . 
The derivation of the full results is under way [14]. 

In [ p^|p6| we have derived approximate expres- 
sions for the 3-loop MS splitting functions which 
are sufficient for reliable NNLO analyses down to 
x ~ 10 -3 . These functions turn out to be much 
less important than the 2-loop coefficient func- 
tions at x > 10~ 2 . Thus it is possible, based on 
partial results [|7p|, p^ , |l8| ] on the 3-loop coefficient 
functions, to proceed to NNNLO at large x pjj , 
especially for the non-singlet case most important 
for extractions of a, from DIS. 
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2. Parton densities: formalism 

It is convenient to work with the flavour non- 
singlet (NS) and singlet (S) combinations of the 
(anti-)quark and gluon densities, qi, qi and g: 



(1) 



Here Nf is the number of effectively massless 
flavours. As in (Q) we often suppress the depen- 
dence on the momentum fraction x and the renor- 
malization and factorization scales, fj, r and (if. 

Using (0) the evolution equations are decom- 
posed into 2Nf — \ scalar (NS) equations and the 
2x2 singlet system, all schematically written as 
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At fi r = fj,f the expansion of the splitting functions 
V up to NNLO is given by 

a 3 s P^(x) . (3) 
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Our choice of the expansion parameter reads 

a s = a s /4n . (4) 

The expression for fx r ^ (if is obtained from (JsJ) 
by inserting the expansion of a s (^) in terms of 
a s (/LtJ?) and Lr = ln(/ij/^). Large logarithms in 
V are avoided by choosing jj, r = O(fif). 
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3. Splitting functions 

The functions and pW in @ are known B . 

(2)3: 

The current information on P^ s (a;) comprises 

(2) + 

• the first five even-integer moments of P^ s 
given by P$$ + (N) = fjdx x N ^P^ + (x) 

and the first moment (iV=l) of P^", 

• the complete iV| contribution 

• the leading small- a; terms oc In 4 a; JToj . 

(2)— (2) + 

The difference of P^g and P^,g is expected to 
be negligible at large x. It has been conjectured 
that the leading large- a: terms are oc 1/[1— a;] + j2(|. 

The following partial results have been derived 

(2) 

so far for the singlet splitting functions P^ (x) : 

• P}p(N) for N = 2, 4, 6 and 8 §, 

• the C A Nj contribution to Pgf{x) Q, 

• the leading small- a; terms oc (1/x) In a; of 
P%\ P$ II and P$ 0], see also @. 

(2) (2) 

The 1/[1— x]+ terms of P g y and P g q are expected 
to be related by a factor Ca/Cf — 9/4. 

We have derived approximate expressions for 

(2)i (2) 

P^g (a;) and P>- ' (x) from these constraints. 
After decomposing the functions into 

p(2) = pW+NfPP+NfpW , (5) 

we employ the ansatz (cf. ||) 

Tim 

P^>{X) = Y. A Mx)+fe{x) . (6) 

Tl = l 

The basis functions /„ are build up of 1/[1 — + , 
<5(1 — a;) and of powers of ln(l — a;), x, and 
In x. The coefficients A n are determined from the 
n m — 5 (n m — 4) linear equations provided by the 
non-singlet (singlet) moments of after tak- 
ing into account the other constraints collected 
in / e in (^). The remaining uncertainties are es- 
timated by 'reasonably' varying the choice of the 
basis functions /„, typically considering some 20 
to 40 trial functions. Finally two approximations 



spanning the error band are selected, except for 
the highest unknown ^-contributions in (^) for 
which one central representative is sufficient. 
This procedure is briefly illustrated in Fig. FA 

(2) + 

for the Nf = part of P^ s (a;) . The upper plot 
shows 24 trial functions. The approximations A 
and B emphasized in the plot have been selected, 
after considering also the convolution with a typ- 
ical input shape shown for these two functions in 
the lower plot. As can be inferred from Fig. [jj 
the presence of the convolution in (|^) consider- 
ably increases the effective accuracy of our ap- 
proximations illustrated in Figs. [I] and ||: The 
convolutions smoothen out the oscillating large- a; 
differences between different approximations to a 
large extent. They also partly compensate the 
large small- a; uncertainties of P^ present despite 
the x — > constraints of |I(||l2|,[l3| . 
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Figure 1 . Top: P ^J for Nf = 0, as obtained from 
the results of by means of (|6|). Bottom: 

Convolutions of the selected approximations A 
and B with a typical non-singlet input shape. 
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Figure 2. The approximations selected for py 
for Nf — 4. The pure singlet (PS) contribution 

Pqf 1 — P^s + is a l so shown at x > 0.3. 



4. Parton densities: results 

We illustrate the impact of the NNLO terms 
on the parton evolution by the derivatives q = 
d\nq/dln q = q-^s> S, g, at a reference scale 

2 



2 2 



30 GeV^ 



(7) 



corresponding to a s (p^ — (jZ ) = 0.2. The input 
densities adopted for the non-singlet case read 



xq m 



x°- 5 (l -xf . (8) 

For the singlet distributions we employ 

xE(x,fj, 2 fi0 ) = 0.6x-°- 3 (l~x) 3 - 5 (l + 5x°- s ) 

xg(x,n 2 fi0 ) = 1.0a;- - s7 (l-a:) 8 . (9) 

The dependence of the results on the renormal- 
ization scale is presented via 

max q — min q 



Aq EE 



(10) 



2 | average q \ 
where fi r is varied over the conventional interval 
1/4 < £ < Atf . (11) 



The NNLO effects on the derivatives q and the 
NLO and NNLO scale uncertainties Aq are shown 
in Figs. |^ and |]. The present inaccuracies of 
the NNLO results caused by the uncertainties re- 
maining for the functions P^ are represented by 
the bands spanned by the NNLOa and NNLOs 
curves. The central results ±(NNLO A + NNLO B ) 
are not shown separately. 

The uncertainties of the NNLO derivatives q 
due to the approximations for P^ are entirely 
negligible for x £ 0.1. They increase towards very 
small values of x, but do not exceed ±2% above 
x ~ 10~ 3 (or a few times this number for scales 
H f much smaller than (0)). Given the small 
size of the NNLO corrections and the weak fi r - 
dependence remaining at NNLO, one can safely 
estimate that contributions beyond NNLO affect 
the parton evolution, for a s ^ 0.2, by less than 
1% at large x and 2% down to x ~ 10 -3 . 
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Figure 3. Top: The NNLO corrections for q^ s = 
d hi q^ s /d In fjlj at \ij — /x/ ; o for the input (||). 
Bottom: The relative /^-uncertainty of the NLO 
and NNLO results for q^ s using (|To|) and (ITT 
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Figure 4. As Figure [|, but for the singlet quark 
and gluon densities E and g given in (||) and (||). 
The spikes close to x = 0.1 in both figures are due 
to zeros of the respective denominators. 



5. Structure functions: formalism 

The unpolarized non-singlet (a = 1, 2, 3) and 
singlet (a = 1, 2) structure functions F a are ob- 
tained by convoluting the solutions g^ 2 ,^ 2 .) 01 
(||) with the corresponding coefficient functions, 

77 a F Q (a;,Q 2 ) = [C a (a s , L M , L R ) <Z> q](x) . (12) 

For L M = In (Q 2 /M/) 7^ 0, but Mr = /< / = 0) : 

(13) 

with C a — ( C a(? , Can ) in the singlet case. The 
coefficients r\ a in (|l3|) include the charge factors 



C a = ci°)(x)+X;ai{c«(x) 



;=i 
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so that c 
„(°) 



(0) 

a,NS ■ 



„(0) 



= 5(l—o;), whereas, of course, 

(Z,m) 



0. The contributions t4' ; fixed by 
renormalization-group constraints are build up of 
the c a k ^ and the splitting functions P^ up to 
k = l — l. The generalization of ( fi"3] ) to /i r ^ /// 
proceeds as indicated below (0). 



The scaling violations of the non-singlet struc- 
ture functions can be conveniently expressed in 
terms of these structure functions themselves 
(thereby removing any dependence on fj,f), viz 

-TTT^Fa,NS = ICa.NS ® Fa.NS ■ (14) 

The kernels /C a ,NS are derived by differentiating 
( p^| ) with respect to InQ 2 and then eliminating 
the quark densities using the same equation. 

6. Coefficient functions 



Besides the functions in (|13[), see 



.,(2) 



also 

are known MM- 



the NNLO contribution c, 
Those expressions are rather lengthy and involve 
higher transcendental functions. We have thus 
provided compact approximations which are suf- 
ficiently accurate for any foreseeable application. 

As illustrated below (Fig. 6), the impact of the 

( 21 

functions c a (especially of the quark coefficient 
functions which contain large soft-gluon emission 
terms) is much larger than that of the splitting 
functions P^ at x > 10~ 2 . The same situation is 
expected for the NNNLO quantities c a ^ and P^ . 
Hence a good approximation to the NNNLO at 

(3) 

large x can be obtained by just retaining the c a ■ 

(3) 

The current information on the c a comprises 

• the first five even-integer moments of c^ 3 ^ , 
and the first four of and c^ g @|§, 

• the four leading large- a: terms oc ln fc (l— x)j 



[1 -*]+, * = 2, ...,5 of eg 
IQ, fixed by the results of |l 
with those of [g for k = 2 

(3)- 



and 



,(3) 



together 



For c 2 only the first moment (= 0) is known 
from the Adler sum rule. However, the results on 

(3) 



(2)± 

c 2 ns indicate that the difference <■ 



-r (3)+ has 
2,NS L 2,NS 



a negligibly small effect. Results for the lowest 
moments of C3 will become available soon |2^| . 

Focusing on the non-singlet case most relevant 
for cn s -determinations from structure functions, 
we have employed the above information to derive 
approximations of Cjjjgfa:) [@- The impact of 
their residual uncertainties is small for x > 0.2. 
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7. Structure functions: results 

We illustrate the effect of the NNLO (and 
NNNLO) terms on the structure functions at 

Q 2 w 30 GeV 2 . (15) 

In ( |l4| ) we employ the non-singlet input shape 

(x,Q 2 ) = x°- 5 (l-x) 3 . (16) 



F + 

r 2,NS 



For the singlet case we fix, besides a s (Q 2 ) = 0.2, 
the input parton densities (^), hence not F 2t s- 
The /^-dependence of the results for / = F 2t s, 
^2,NS = dhii^Ns/cHnQ 2 , F 2 ,s and F 2 ' s = 
dF 2: s I 'd\nQ 2 is presented via 

. ■ max f — min f 1 , , . . 

jf, tQ 2 <4Q 2 . (17) 

2 | average / | 4 



The singlet results are shown for [if = \i T 
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Figure 5. Top: The NNLO and NNNLO cor- 
rections for F 2 + NS for the input ( |l5| ) and (|l6|). 
Bottom: The relative /^-uncertainties ( |l7j ) com- 
pared to the corresponding NLO result. 



The results for p2,NS are shown in Fig. ||. The 
uncertainty bands for the NNNLO predictions 
take into account both the remaining inaccura- 



cies of the coefficient functions c 2 n S (x) and the 
possible effects of the splitting functions P^ s . At 
0.25 < x < 0.7 the // r -uncertainties Ai<2,NS are re- 
duced by a factor of two (four) or more at NNLO 
(NNNLO). These uncertainties lead to the follow- 
ing estimates for the errors of a s (M§) due to the 
truncation of the perturbation series: 



Aa s (M|) NLO 
Aa s (M 2 ] 
Aa s (M|; 



NNLO 



NNNLO 



0.005 
0.002 
0.001 



(18) 



A 1% accuracy is achieved at the NNNLO level. 
As the scaling violations for the same a s (Q 2 ) are 
stronger at NNLO and NNNLO than at NLO, 
higher-order fits of data on i^NS wu l yield some- 
what lower central values of a s (M§ ), 



a s(A^f)NNNLO — a s (M§)NLO 



-0.002 



(19) 



The results for the singlet case are presented 
in Fig. ^. F 2 ,s receives large positive corrections 
at large x, caused by the soft-gluon parts of the 
quark coefficient functions. The sizeable negative 
NNLO corrections at small x are dominated by 
the gluon contribution. It is worth noting that 
the positive l/x term of Sp g does not dominate 
this correction even at x < 10~ 3 . 

The Q 2 -derivative of F 2j s is dominated by the 
quark contribution at x > 0.3, and by the gluon 
contribution at x < 0.03. Thus we present the 
logarithmic derivative i*2,s in the former x-range, 
and the linear derivative F 2 s in the latter region. 
Note that the positive NNLO gluon contribution 
reaches 5% of the total \F 2 .s\ at x = 0.5, enough 
to jeopardize purely non-singlet analyses of F 2 P 
data also in the region x > 0.3. 

The reduced /^-dependence of both F 2 and its 
derivatives leads to a better theoretical accuracy 
of determinations of the parton densities from 
data on F 2 ,s and dF 2 ,s/d\nQ 2 at Q 2 ~ 30 GeV 2 : 
NNLO uncertainties of less than 2% from the 
truncation of the perturbation series are obtained 
for the quark density at 10~ 3 < x < 0.5 and for 
the gluon density at 3 • 10 _3 < x< 0.2. 
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Figure 6. Top: The NNLO corrections for F2 = 
F2.S and its Q 2 derivatives for the input (||) at 
A 4 / = Q 2 ■ Bottom: The corresponding (i r = (if 
uncertainties (|l^) compared to the NLO results. 

8. Summary 

We have briefly discussed the evolution of un- 
polarized parton densities and structure functions 
in the MS scheme. Our approximate results for 
the 3-loop splitting functions P^ 2 \x) pave the 
way for promoting, even though only at x> 10~ 3 , 
global analyses of DIS and related processes to 
NNLO accuracy. We will also provide approxima- 
tions for the 3-loop non-singlet coefficient func- 

(3) 

tions c^g, thus enabling NNNLO determinations 
of a s from structure functions at least at x>0.2. 

At very large x, x > 0.8, terms even beyond 
NNNLO are relevant. Here results are available 
from soft-gluon resummation [p7|Jl8[ . Progress 
towards the important HERA small- a; region of 
x < 10 -3 at moderate/low Q 2 requires the full cal- 
culation of the three-loop splitting functions Q] . 

Fortran subroutines of our parametrizations 
of the 2-loop coefficient functions and our ap- 
proximations of the 3-loop splitting functions can 
be obtained from neerven@lorentz.lcidcnuniv.nl 
or avogt@lorentz.leidenuniv.nl. 
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